Weak subordination breaking for the quenched trap model 
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We map the problem of diffusion in the quenched trap model onto a new stochastic process: 
Brownian motion which is terminated at the coverage "time" S a = '^2°^-_ 00 ( n x) a with n x being 
the number of visits to site x. Here < a = T/T g < I is a measure of the disorder in the original 
model. This mapping allows us to treat the intricate correlations in the underlying random walk in 
the random environment. The operational "time" S a is changed to laboratory time t with a Levy 
time transformation. Investigation of Brownian motion stopped at "time" S a yields the diffusion 
front of the quenched trap model which is favorably compared with numerical simulations. In the 
zero temperature limit of a — > we recover the renormalization group solution obtained by C. 
Monthus. Our theory surmounts critical slowing down which is found when a — > 1. Above the 
critical dimension two mapping the problem to a continuous time random walk becomes feasible 
though still not trivial. 

PACS numbers: 05.40. Jc,02.50.-r,05.20.-y, 46.65. +g 
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I. INTRODUCTION 



Random walks in disordered systems with a diverg- 
ing expected waiting time have attracted vast interest 
over many decades. Two approaches in this field are the 
annealed continuous time random walk (CTRW) model 
and the quenched trap model (QTM). Starting in the 
70's, the Scher-Montroll CTRW approach was used to 
model sub-diffusive photo-currents in amorphous mate- 
rials [H-[3] and for contaminants transport in hydrology 
Q . Bouchaud showed that the trap model is a useful tool 
for the description of aging phenomena in glasses 
Then fractional kinetic equations which describe CTRW 
dynamics became a popular tool [9( . More recently these 
models were used to describe non self averaging [H| 
and weak ergodicity breaking [f| [l2| which is important 
for the statistical description of dynamics of sing le q uan- 
tum dots [lH and single molecules in living cells [3 . 

This manuscript presents a new approach for random 
walks in a quenched random environment i.e. site disor- 
der at each lattice point is fixed in time. In its generality 
this topic has attracted tremendous interest in Physics 
ElHil and Mathematics [i^Sl. For the QTM the 
critical dimension is two [l?], l26l - [29| . Above two dimen- 
sions the Scher-Montroll continuous time random walk 
(CTRW), which is a mean field theory, qualitatively de- 
scribes the sub-diffusive process. According to Polya's 
theorem |30L l3l| on a simple lattice and in dimension 
three, a random walk is non recurrent. Hence in a dis- 
ordered system the particle (roughly speaking) tends to 
visit new lattice points along its path. In contrast in one 
dimension the random walk is recurrent, and a particle 
visits the same lattice point many times. Thus above the 
critical dimension the CTRW approach works well, but 
fails in one dimension, due to correlations of the ran- 
dom walk with the disorder. In other words renewal 
theory used within the annealed CTRW framework is 
not a valid description of the QTM [17] . Beyond mean 
field renormalization group methods are used to tackle 



the problem of random walks in quenched environments 
[HIM [H HI]. For example Machta Hlfound the scal- 
ing exponents of the QTM and Monthus [32[ investigated 
its diffusion front in the limit of zero temperature (see de- 
tails below). While the renormalization group method is 
powerful, it has its limitations: a simple approach which 
can predict the diffusion front of random walkers in the 
QTM is still missing. 

We provide a new approach for random walks in the 
QTM which we call weak subordination breaking. For 
CTRW it is well known that one may decompose the 
process into ordinary Brownian motion and a Levy time 
process, an approach called subordination |34Ho8l [. In 
this scheme normal Brownian motion takes place in an 
operational time s. The disorder is effectively described 
by a Levy time transformation from operational time s 
to laboratory time t (see details below). This method 
is not intended for random walks in fixed random envi- 
ronments since it is based on the renewal assumption. 
The latter implies the neglect of correlations in the sense 
that waiting times are not specific to a lattice site. So a 
new approach capable of dealing with quenched disorder 
is now investigated. A brief summary of our results was 
published in [39(. 

This manuscript is organized as follows. After present- 
ing the QTM in Sec. (JTTJ) we briefly review the standard 
subordination scheme in Sec. (|III[) . The concept of ran- 
dom time in the QTM is presented in IIVI which leads 
to weak sub-ordination breaking in Sec [V] General prop- 
erties of the diffusion front (P(x,t)) are found in Sec. 
rVTl while Sec. |VH] and [VTTTI deal with the limits strong 
and weak disorder respectively. Sec. 1 1X1 discusses criti- 
cal slowing down. All along the work we compare theory 
with numerical simulations. 



II. QUENCHED TRAP MODEL [HI, H3, S3, S3 

We consider a random walk on a one dimensional lat- 
tice with lattice spacing equal one. For each lattice site 
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x there is a quenched random variable r x which is the 
waiting time between jump events for a particle situated 
on x. After time t x has elapsed the particle jumps to 
one of its two nearest neighbors with equal probability. 
The particle starts on the origin x — at time t = 0, 
waits for time To, then jumps (with probability 1/2) to 
x = 1, waits there for t\ etc. Note that if the particle 
returns to x = it will wait there again for a time in- 
terval To- The {t x }s are positive independent identically 
distributed random variables with a common PDF ip{T x ). 
The goal of this paper is to find the long time behavior 
of (P(x,t)), the probability of finding the particle on x 
at time t averaged over the disorder. 

In the literature two related models are usually con- 
sidered. The first model, which we use in simulations 
presented below, assumes that for a given lattice site x a 
particle will wait for a fixed waiting time t x . A slightly 
more physical approach is to assume that waiting times 
on lattice point x are exponentially distributed with a 
mean t x . Bertin and Bouchaud [4fJ, showed that the 
two approaches yield the same asymptotic results in the 
limit of long measurement times. 

In this manuscript our main interest is with power law 
waiting times 



-a 



-(l+a) 



(1) 



for t — > oo and < a < 1. The mean waiting time (t x ) 
x and in this sense the diffusion is scale free. According 
to Tauberian theorem [30| the Laplace transform of the 
waiting time PDF is 



4>{u) ~ 1 - Au a + 



(2) 



when u — > 0. In the QTM the physical mechanism lead- 
ing to these power laws is based on trapping dynamics 
[17| . On a lattice points a; we randomly assign traps. The 
energy depth of the trap on x is E x and the process of 
activation from a trap is thermal. According to Arrhc- 
nius law t x cx exp(E x /T) where T is the temperature. 
Then assume that the PDF of E x > is exponential 
f(E x ) = ex.p(—E x /Tg)/T g where T g is a measure of the 
energy disorder. One easily finds 



T 



and A = |r(-a)|a. 



(3) 



Due to the Boltzmann factor t x cx exp(_E/T) small 
changes in energy lead to exponential changes in waiting 
times, thus it is enough to have an exponential distribu- 
tion of energy traps to obtain power law waiting times. 
Experimental observation of the linear dependence of 
a on temperature, in photocurrent spectroscopy in a - 
As2Se3 can be found in Fig. 3 in 42j. We note that the 



stochastic dynamics under investigation describes sev- 
eral other mechanisms of anomalous diffusion beyond the 
QTM [l3]. For example random walks on comb struc- 
tures with power law distributed lengths of the comb's 



teeth which mimics a random walk on the percolation 
cluster. Thus the trap model describes both energetic 
disorder and spatial disorder. 

In what follows we will also consider the limit a — > 
1. This limit is meant in the sense that ip{u) ~ 1 — 
Au.... which means that the average waiting time is finite 
(Gaussian diffusion front). The very special border case 



■>P(t x 



oc t was treated by Bertin and Bouchaud 4C 



It yields Gaussian diffusion with logarithmic corrections 
and is not treated here. 



III. SUBORDINATION IN THE ANNEALED 
TRAP MODEL (=CTRW) 

We now briefly review the annealed version of the 
model: the well investigated Scher-Montroll- Weiss con- 
tinuous time random walk (CTRW) 0, [13, H3| in partic- 
ular we discuss the concept of time subordination [34H38I] . 
Later we contrast the CTRW approach with the intricate 
problem of the quenched type. The CTRW model con- 
sidered here is for a one dimensional random walk on a 
lattice with lattice spacing equal unity. Starting on the 
origin x = at time t = the particle waits for a time 
t\ it then jumps to one of its nearest neighbors (lattice 
points x — +1 ot x — — 1) with equal probability. The 
process is then renewed, namely the particle waits on 
lattice point +1 (for example) for time t% until it jumps 
back to or +2 etc. The waiting times {t\,t2, • • • , t n , ■ ■ •} 
are independent, identically distributed random variables 
with a common PDF <j) a (t). Here t n is the nth waiting 
time, which is not correlated with a specific lattice point 
x and hence clearly the CTRW model is very different 
from the quenched case. Similar to the quenched case we 
consider waiting time PDFs with a diverging averaged 
waiting times / °° t(j) a (t)dt = oo namely 



K (t) 



A, 



\n- a )\ 



-(l+a) 



(4) 



with < a < 1 and A a > 0. The corresponding Laplace 
transform of the waiting time PDF behaves like 



4> a {u) - 1 - A a u a + 



(5) 



when u is small. As well known the diffusion is anomalous 

(x 2 ) (xt a mini. 

Let P(x, t) be the probability of finding the particle on 
x at time t. By conditioning on the number of jumps s 
performed till time t 



p (x,t) = ^2n t (s)q(x,s) 



(6) 



s=0 



where n t (s) is the probability of performing s jumps in 
time interval (0, t) and q{x, s) is the probability that after 
s steps the particle is located on x. In the limit of large 
t the number of jumps s is also large. Following 17] we 
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apply the Gaussian central limit theorem 

1 



q(x,s) 



-e 2s 



(7) 



which is valid when s — > oo. Here we used the model 
assumption that the variance of the jump lengths is unity 
i.e. the lattice spacing is equal 1. 

To find nt(s) we consider the random time 



(8) 



In the limit of large s the time is a sum of many inde- 
pendent identically distributed random variables with a 
diverging mean waiting time (since < a < 1). Hence 
Levy's limit theorem applies. Let 



t 



V 



,1/q 



(9) 



then in the s — > oo limit 



= 1- 



A, 



(10) 

Namely the PDF of rj > is a one sided Levy function 
denoted with l a ,A a ,i.( r l) w hich is defined via its Laplace 
pair 



-Aa 



(11) 



These Levy PDFs are well investigated: their series ex- 
pansion, asymptotic behaviors, and graphical presenta- 
tions can be found in RTi 1 ■") . Information on these PDFs 
essential for our work are summarized in Appendix A. 
From the PDF of 77 we find the PDF of s. Since both s 
and t are increasing along the process the transformation 
is straight forward. Using Eq. (ITU)) and ?/~ Q = s/t a we 
find the well known PDF of s |l7j 



n t (s) 



t 



t 



,1/a 



(12) 



In the long time limit the Green function of the CTRW 
process is thus given by [13, l34j 



P(x,t) 



n t (s)- 



-.ds 



(13) 



where we switched from a summation in Eq. to inte- 
gration. 

The time transformation Eq. (|13p maps normal Gaus- 
sian diffusion to anomalous diffusion. In [35j j P(x, t) was 
obtained in d dimension by solving the integral transfor- 
mation which applies more generally to solutions of the 
fractional time Fokker-Planck equation (46[. More im- 
portantly, we may think of s as an operational time in 
which usual Brownian motion is performed. The oper- 
ational time s is a random variable whose statistics is 



determined by the PDF n t (s) where f is a laboratory 
time. In other words the annealed disorder turns the 
operational time to a random variable. We note that 
subordination scheme can be formulated for the trajec- 
tories of the corresponding paths, in a continuum limit of 
the walk and was the topic of extensive research (47l - [54j . 

Not surprisingly subordination of this type does not 
work for the QTM in one dimension. As mentioned 
in the introduction the process of a random walk in a 
QTM is clearly not a simple renewal process. The parti- 
cle returning to a lattice point visited already, "remem- 
bers" its waiting time there. Mathematicians have rig- 
orously shown that in dimensions higher than one [281 ] 
or in the presence of a bias [55[ (see also [56U58j ) the 
CTRW approach describes well the quenched dynamics 
since the particle does not tend to revisit the same lat- 
tice points many times, thus confirming physical insight 
in [Tjl [26[ (for dimension d — 2 logarithmic correc- 
tions are also important). While the three dimensional 
QTM belongs to the domain of attraction of the CTRW 
the calculation of the anomalous diffusion constant is not 
trivial (see discussion in the summary). Here we focus 
our attention on the unsolved case: the QTM in one di- 
mension since there the Scher-Montroll CTRW picture 
d G3, H| breaks down. 



IV. TIME IN THE QUENCHED TRAP MODEL 



The time t in the QTM is 

OO 



(14) 



where n x is the number of visits to lattice point x which 
we call the visitation number of site x. Since we are inter- 
ested in (P(x, t)} where the brackets are for an average 
over the disorder, we will consider ensembles of paths on 
a large ensemble of realizations of disorder. As mentioned 
the {r^js are independent identically distributed random 
variables with a common PDF iI)(t x ) and the {n x }s are 
also random variables. 

Let us consider the random variable 



t 



where 



(15) 



(16) 



and we call S a the a coverage time. At this stage it is 
convenient to consider paths where S a is fixed and t is 
random and later we will switch to the opposite situa- 
tion (similar to the arguments for s and t in the CTRW 
model). When a — 1, S a is the total number of jumps 



made J2 X = 



s. In the opposite limit a — > 0, the 
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a coverage time Sq is the distinct number of sites vis- 
ited by the random walker which is called the span of the 
random walk. Notice that t in Eq. (fl"4|) is a sum of non 
independent and non identical random variables. 

We show that the PDF of 77, in the limit S a — > 00 is a 
one sided Levy stable function 



The PDF of i] is: l a<AA (77) . 



(17) 



Namely the heavy tailed distribution of the waiting times 
t x determines the statistics of 77 through the characteris- 
tic exponent a, while the visitation numbers {n x } deter- 
mine the scaling through S a . By definition the Laplace 
rj — > u transform of the PDF of 77 is 



(e-""> = (exp 



E'Hn 
TJ-'< 

t=-oo 



(18) 



We average with respect to the disorder, namely with 
respect to the independent and identically distributed 
random waiting times t x , and obtain 



-um 



n° 



(s a y/<* 



(19) 



where ip(u) is the Laplace transform of the PDF of wait- 
ing times iI){t x ). Now assume ip( u ) = ex P( — Au a ) ~ 
1 — Au a H . Then using Eq. (|T9|) we have 



IE 



,exp 



A(n x ) a u a 



S a 



(20) 



Hence if the waiting PDF is a one sided Levy PDF, i.e. 
ijj(u) = exp(— Au a ), so is the PDF of r\. In Appendix B 
we consider the general case where ip(r x ) belongs to the 
domain of attraction Levy PDFs (i.e. families of PDFs 
satisfying Eq. ip(u) <~ 1 — Au a + •••). We there prove 
that the statement in Eq. (flTf is valid. 

From Eq. (fTT]) we learn that the CTRW operational 
time s, that is the number of jumps made in the random 
walk, looses its importance in the quenched model. In 
the QTM the operational time is the a coverage time S a . 
We now invert the process fixing time t to find the PDF 
of S a 



n, 



(S a ) ~ — (S a ) 1//a 1 l a .A.i 
a 



(21) 



In the next section we explain how to use the operational 
time S a to obtain the desired diffusion front of the QTM. 



V. WEAK SUBORDINATION BREAKING 

To find the solution of the problem, namely find 
(P(x,t)) for the QTM we follow the following steps: 

1. Choose the laboratory time t which is a fixed pa- 
rameter. 



2. Use a random number generator and draw the sta- 
ble random variable rj from the one sided Levy law 
la,A,i(n) . 

3. With 77 and t determine the hitting target S a which 
according to Eq. (fT5)) is S a = (t/n) a . 

4. Generate a Binomial random walk on a lattice, with 
probability 1/2 for jumping left and right. Stop the 
process once its S a crosses the hitting target set in 
step 3. 

5. Record the position x of the particle at the end of 
the previous step. 

6. Go to step 2. After this loop is repeated many 
times, we generate a histogram of x. 

The histogram once normalized yields (P(x, t)} when t 
is large. Notice that in this scheme there is no disorder. 
The second step is implemented with a simple algorithm 
provided by Chambers et al These authors show 

how to generate stable random variables like 77 using two 
independent uniformly distributed random variables and 
for convenience their formula is provided in Appendix A. 

More importantly we can now start treating the prob- 
lem analytically, and find the diffusion front. So far we 
have replaced the problem of random walks in the QTM, 
with a new stochastic process: Brownian motion which 
is stopped when the hitting target S a is crossed. In other 
words we got rid of the disorder. Notice that the CTRW 
process and standard subordination j3^ - l38j are reached 
once we replace S a with s. In this sense the QTM ex- 
hibits what we call weak sub-ordination breaking: the 
operational time is now S a still the Levy time transfor- 
mation used already in the usual subordination scheme 
Eq. (fTB"]) remains a useful tool. 



VI. THE DIFFUSION FRONT OF THE 
QUENCHED TRAP MODEL 

Let Pb(x, S a ) be the PDF of x for a binomial random 
walk on a lattice at operational time S a . The subscript 
B indicates that the underlying motion is Brownian. The 
corresponding paths are generated from a random walk 
on a one dimensional lattice, with equal probability of 
jumping left and right, which is stopped when S a is 
reached (or crossed for the first time). The averaged over 
disorder propagator of the QTM is found using Eq. (|2"Tj) 
and the scheme presented in the last section: 



(P(x,t)) 



P B (x,S a )n t (S a )dS a , (22) 



which is valid in the long time limit. Eq. (|2"2"|) is a gen- 
eralization of the subordination equation (|13|) . Namely 
it transforms Brownian motion stopped at the coverage 
time S a to the QTM dynamics in laboratory time t. 

From Eq. (|2"21 we may find general properties of the 
Green function (P(x, t)) in terms of its corresponding 
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FIG. 1: The PDF B a (z) exhibits a transition between a 
Gaussian shape when a — > 1 to a V shape when a — > 0. 
Simulations of Brownian motion on a lattice yield excellent 
agreement with theoretical predictions Eqs. (|39l I50|) without 
fitting. 



Brownian partner Ps(x,S a ). For example the Laplace 
t — > u transform 



(P(x,u)) = Au^PBix^u ). 



(23) 



Less formal relations are found if we exploit the scaling 
behavior of Brownian motion as now explained. 



A. Scaling arguments 



Brownian motion follows the usual diffusive scaling 
x 2 oc s where s is the number of steps. In Appendix C we 
show that s cx (<S Q ) 2/(1+Q ° which is now explained using 
simple arguments. For Brownian motion the particle ex- 
plores a region which scales like s 1 / 2 . The visitation num- 
ber n x within this region (i.e. roughly \x\ < s 1 / 2 ) is the 
number of jumps made (s) divided by the number of sites 
in the explored region (s 1 / 2 ) so n x oc s/ s 1 / 2 = s 1 / 2 . Since 
particles typically visit \x\ >> s 1 / 2 only rarely n x oc 
there. Hence S a cx \fs(n x ) a oc s^ 1+a ^ 2 . Indeed in Ap- 
pendix C we show that 



(Sa) — C a 



with 



r — 



(24) 



(25) 



When a = 1 we have Ci = 1 since 5i = y~!°° n T = s. 
In the opposite limit a = we find a well known result 
obtained by Dvoretzky and Erdos [6(3] 



(So) 



(26) 
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B a (z = 0) 
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0.1 


0.592 


0.08 


0.2 


0.673 


0.15 


0.3 


0.746 


0.2 


0.4 


0.808 


0.24 


o..-) 


0.859 


0.28 


0.6 


0.907 


0.3 


0.7 


0.929 


0.33 


0.8 


0.961 


0.35 


0.9 


0.986 


0.38 


1 


1 


1/V27T 



TABLE I: Brownian simulations on a lattice give (z 2 ) and 
B a (z — 0) which in turn provide the corresponding solution 
of the QTM with Eqs. (|32l 133) . 



By definition (Sq) is the averaged number of distinct sites 
visited by an unbiased random walker [30j . 

Using x oc s 1 / 2 and S a oc s( 1+a )/ 2 scalings we have 
x oc (S a ) 1 ^ 1+a K We emphasize that this is a property of 
simple binomial random walks which we can now exploit 
to investigate the solution of the QTM. More specifically 
this scaling implies 



P B (x,S a ) = 



(5 Q ) 1 /(i+«) 



B a 



( l S Q ) 1 /(i+«) 



(27) 



Here B a (z) is a a non negative function normalized ac- 
cording to J_ oo B a (z)dz = 1. Further from symmetry 
of the walk B a (z) — B a (—z). As shown in Fig. [TJ the 
PDF B a (z) exhibits an interesting transition between a 
V shape for a — > and a Gaussian shape when a — > 1. In 
the following sections we will investigate B a (z) in detail. 

With B a (z) we obtain useful relations between the dif- 
fusion front of the trap model and Brownian motion. De- 
fine the dimensionless time t = t/A l / a and the scaling 
variable 



(28) 



(t)n 

Then it is easy to show that 



(P{*,t)) ~ %Q 



(29) 



and using Eq. fl22J) 

POO 

g a (0= / dyy^B a ^y^)l aX1 (y). (30) 
Jo 

For the behavior of (P(x, t)) on the origin we use 

if q/a > 1 

if q/a < 1 



r f 00 

/ dy y n aAA (y) = { 
Jo rq-g/a) 
I r(i- Q ) 



(31) 
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QTMd=0.4 


<> 


QTM 0=0.6 




QTMd=0.8 




Eq. (33) 




Eq. (53) 




Eq. (47) 




FIG. 2: The mean square displacement versus time. Numer- 
ical data obtained from direct simulations of the QTM match 
perfectly the theory based on weak subordination breaking 
(the lines plotted using Eq. (|33[l and Table [XJ . Analytical 
predictions for (a; 2 ) Eq. (TJTJ) for a = 0.2 and Eq. for 
a — 0.4, 0.6, 0.8 perfectly match direct simulations and weak 
subordination scheme (based on Brownian simulations) and 
hence the corresponding curves fall one on top of the other. 



and then find 



(P(x = 0,t))~B a (0)- 



(it a) 



(32) 



This is a useful result since the behavior of B a (z) on the 
origin z — gives the corresponding behavior of (P(x = 
0, £)) without the need to solve any integral. Further Eq. 
(I3"2l hints to an interesting behavior when a 0. The 
ratio of the T functions diverges in that limit, hence as 
shown in Fig. [T]i? ct (z = 0) must go to zero when a — > 
for (P(x = 0,t)} to remain finite. Such a behavior is 
analytically investigated in the following section. 

Another useful relation is found between the moments 
(|x|«) = (J^ \ x \iP(x,t)dx) of the original QTM and the 

moments = \z\ g B a (z)dz, Using Eqs. (|30|[3T|) 

we find 



(\x\ q ) = (\z\ 



(33) 



The scaling x 2 oc (J) 1+ ° was obtained long ago in [l7ll27| 
using elegant scaling arguments and in [261 ] using renor- 
malization group approach. The new content of Eqs. (1301 
I32I33P is that once we obtain B a {z) either from theory 



or simulations of Brownian trajectories, we have a use- 
ful method to obtain exact statistical properties of the 
diffusion front. 

On a computer our approach is very useful. For exam- 
ple we have numerically generated Brownian trajectories 
on a lattice and obtained B a (z) in Fig. [T] while (z 2 ) and 
B a (0) are reported in Table Q] With (z 2 ) given in Table U 
and Eq. (|3"3")l we get the mean square displacement of the 
QTM (x 2 ). Direct simulations of the QTM are favorably 
compared with the predictions of our theory in Fig. [2j 

Finally the cumulative distribution function G a (£ < 

S) = /"„ ffa(£)df, the probability that the random vari- 
able £ attains a value less then H is found using Eq. (|2"T)l 



poo 

G a (£<|3|)=l - / dzB a (z)L a 
Jo 



(34) 



Here L a (y) = Jq l a ,i 1 i(y)dy is the cumulative distribu- 
tion of a one sided stable random variable. From sym- 
metry G„(£,-|H|) = 1 - G a (£ < |S|). The integral 
representation of the distribution L a (y) can be found in 



VII. LIMIT a 







As mentioned in the introduction the diffusion front 
(P(x,t)) was treated using a renormalization group 
method by Monthus [32j]. We will now investigate this 
interesting limit using our new approach. For that we 
must first find B a (z) in the limit a — > 0. 



A. lim a _>o B a (z) has a V shape 

We consider Brownian motion on a lattice which is 
stopped once the distinct number of sites visited by the 
walker reaches the threshold Sq and as a reminder Sq 
is called the span. The position of the particle is then 
x and we are interested in the probability P(x,Sq) of 
finding the particle on x. 

The particle starts on the origin, hence clearly we have 
\x\ < Sq. Further P(x — 0, <5>o) = 0, since a particle 
starting on the origin cannot reach the threshold Sq when 
it is on the origin: i.e. a particle returning to the origin is 
not increasing Sq since the origin is not a new site visited 
by the walker. From symmetry P(— x,So) = P(x,So). 

Consider first P(x = So, So). After the first step the 
particle can be either on x = 1 with probability 1/2 
or on a; = —1 with the same probability. Clearly a 
trajectory going through x — —1 cannot contribute to 
P(x = So, So) since to reach x — Sq through x = — 1 the 
span must be at-least of length So + 1. So we consider 
only trajectories going through x = 1. Trajectories go- 
ing through x = 1 are divided into three non-intersecting 
categories. Trajectories that (i) never reach the origin 
x = along their path, (ii) trajectories that reach the 
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On the other hand if it jumps to x = 1 the distance the 
particle must travel is 5o — 2 and the span is So. Hence 
we have 

P(x = So - 1,5b) = \p(* = So,S )+^P{x =So- 2, S ) 

(36) 

where the first (second) term on the right had side de- 
scribes trajectories starting on the origin and in the first 
step jumping to x = — 1 (x = 1) respectively. Similarly 
for Sq — n > with n > being an integer we have 

P(x = S - n,So) = 



s 

S r -1 



1 



m o 
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o o o 
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FIG. 3: (a) Trajectory with a span So. The particle returns 
to the origin along its path [category (ii) in text] , hence in the 
time interval s > 1 (excluding the first step) the displacement 
is So — 1 and the span So. (b) A random walker with span 
So reaching x = So — 1 must pass through x = — 1. Here the 
first jump event brings the particle to x = — 1 and hence the 
span in s > 1 is So and the displacement So- 



origin but never cross it see Fig. 03a) and (iii) trajecto- 
ries that go below x = 0. The latter will have a total 
span greater than So and hence do not contribute. For 
class (i) the span (after stepping into x = 1 in the first 
step) is So — 1. Similarly for class (ii) the span is So- For 
both cases the displacement (from x = 1 to So) is clearly 
So — 1. Hence we have 



P(x = S ,S ) 



~[P(x = S - l,S ) + P(x =S -1,S - 1)] (35) 

where the first (second) term on the right hand side de- 
scribes trajectories returning (never returning) to the ori- 
gin. The half in front of the square brackets is due to the 
displacement in the first jump event. 

Continuing with similar reasoning consider 
P(x = So — l,«So). The particle after the first step 
can be either on x = 1 or x = — 1. As shown in Fig. 
|3Jb), if it is on x = — 1 it must travel a distance So to 
reach its destination So — 1 while keeping the span So- 



^P(x = S -n- l,S ) + ^P(x = S -n+ l,S ) 



Eqs. (|35I37[) are easily solved 



(37) 



P(x,S ) 



So (So + 1) 



for - S < x < S (38) 



and x £ Z. In the limit So 3> 1 we have for the scaled 
variable z — x/Sq the PDF 



lim B a (z) = 



forlzl < 1 



otherwise 



(39) 



We see that B a= o(z) has a V shape. This reflects the 
observation that the particle reaching a large span iSo is 
most likely far from the origin, and the probability of 
reaching the span So for the first time, while the particle 
is on the origin being zero. We now use this property 
of Brownian motion to solve the quenched trap model in 
the limit a — > 0. 



B. Diffusion front in the a — > limit 

Define the Fourier transform of the scaling function 
?«(£) Eq. 



/oo 
-oo 



(40) 



which as usual is also the moment generating function 



q=0 



(2g)! 



(41) 



According to our theory the moments ((£) 2<? ) — 
JZo ff"(C)C 29 dC and similarly (x 2q ) for the QTM are de- 
termined by Brownian motion with the help of the PDF 
B a (z) Eq. (|33|) . In the limit a —> we find using Eq. 
(133 



r 1 1 

? 2q ) = 2 / z 2q zdz = 

1 Jo 1 + 



(42) 
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and hence for a — > Eqs. (1281331) give 

(2?)! 



(a) a=0. 1 



2<?\ _ 



Therefore 



Hm 9Q (^)=£(-ir(-l T ) (fc ? ) 2 «. 



9=0 

summing this series we find 



lim g a (k{) 

a— ^0 



In [1 + (fc 5 ) 2 



Inverse Fourier transform yields 



lim o a (0 = e 

a— 



KIMKI) 



(43) 



(44) 



(45) 



(46) 



where £i(£) = fe°( e t /t)dt is the tabulated exponen- 
tial integral [6l|. This result (written in a different but 
equivalent form) was obtained by C. Monthus 32] using 
the renormalization group method, which is exact in the 
limit a — ¥ 0. In Fig. Hwe show <7 a (£) for simulations 
of the QTM (a = 0.1), Brownian simulations using weak 
subordination breaking outlined is Sec. IVland analytical 
curve Eq. (|4"6"|) . We see that the theory which is exact 
when a —> works well also for small values of a. 

When a is small we find a useful approximation for the 
moments. Inserting (\z\ q ) Eq. (1421) in Eq. (|33|) we have 



(M 9 > 



r ( -2- 

1 l+a 



9 a r 



l+a 



t 



A l/a 



(47) 



Notice that in this limit T[q/(l + a)]/{aT[aq/(l + a)]} ~ 
r(l + q) hence for q = we have (|x|°) = 1 as expected 
from normalization. In Fig. [2] we show (x 2 ) versus time 
for a = 0.2. Numerical simulation of the QTM perfectly 
match Eq. (|4"T|). Note that the theory based on Table 
U and Eq. (f3"3"f does a slightly better job since that ap- 
proach is not limited to the a << 1 regime. 



VIII. APPROACHING THE GAUSSIAN LIMIT 

a = 1 

In this section we consider the case a — > 1 from below. 
We now find an approximation for B a {z) which yields 
the solution of the QTM in this limit. 



A. B a (z) is Gaussian when a — > 1 

As before to find B a (z) we consider Brownian motion. 
The probability of finding the particle on x at time s is 
a Gaussian 



P(x,s) 



exp 



K) 



(48) 




FIG. 4: The diffusion front of QTM (simulations) nicely 
matches our theory based on weak sub-ordination breaking 
[the algorithm in Sec. (|V|) ] and analytical predictions Eqs. 
(HI EU) for a = 0.1,0.75 respectively. 



as is well known. For a — - 1 we have 5i = V!? 3 „ n T = 

1 ^ — ^a;= — oo 

s, namely Si is not a random variable at all since it is 
equal to the number of steps made in the underlying ran- 
dom walk. In other words the PDF of Si is a delta func- 
tion centered on s. Therefore when a is close enough to 
1 we may neglect fluctuations. This means that we omit 
the average in Eq. (1231) and use S a — C a s~^~ . This ap- 
proach together with Eq. ([48]) gives the PDF of finding 
the particle on x for a random walk stopped at the a 
coverage time S a 



exp 



P B (x,S a ) 



2tt (S a /C a ) 



1/2 ' 



(49) 



Hence 



exp 



B a {z) 



and it follows that 



(C a )~z 2 
2 



27T/(C a )^ 



(C a ) 



1/2 



(50) 



(51) 



In Fig. [DBq^z) obtained from Brownian simulations is 
compared with the analytical prediction Eq. ([5TJf for 
a = 0.9. 
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B. (P(x,t)) when a ~ 1 



From Eqs. (|22I49[) we have 
(P(x,t))c 



oo exp 



2{S a /C a ) 1+ 



II 



2tt (S a /C a ) — 



-L n (S a ,t). (52) 



From Eq. (|52|) it is easy to get g a (£,) which is given in 
Eq. fliTTl) in Appendix D. In Fig. @]we compare between 
the scaling function obtained analytically and numerical 
simulations of the QTM, and with Brownian simulations 
according to the disorder free algorithm in Sec. [V] The 
theory works reasonably well even for a — 0.75. 

Using Eqs. (|33I51[) we find the mean square displace- 
ment of the QTM 



<^> ~ (c a y- 



_2_ 1 I !+" J 



t 



(53) 



Calculation of other moments is as simple, since the 
reader may easily obtain (\z\ q ) from the Gaussian PDF 
Eq. flSUJl and then apply Eq. (|33|). The behavior on the 
origin is found using Eqs. (|32|50|) 

(P(, = ll,t)) ^1s (f)-Tfe. (54) 



When a = 1 we get the expected behavior (P(x = 
0, t)) — (27rt) -1 / 2 , which is normal diffusion. 

The scaling function g a (£) is analyzed in Appendix 
D. Using properties of stable PDFs, we show that when 
£«1 



with 

,9a (0) - 



1 + a 



r 

( l-a\ ? 



(C. 



a )^ r (lf^) 



(55) 



(56) 



In the limit a — > 1 we use \im a ^iC a = 1 and Eq. ([55 
gives 



a—>l 



l 



lim 



(57) 



The first term clearly reflects an ordinary Gaussian dif- 
fusion front. The second term vanishes in the limit a = 1 
since T(0) = oo. This is because <7i(£) is Gaussian and 
hence the second term in the expansion must be a £ 2 
term. So the l/r(0) kills the £ Q in Eq. ([57]) as a -> 1. In 
the opposite limit of £ >> 1 a steepest descent method 
gives 



-b 2 C 



(58) 




FIG. 5: We show (z 2 ) versus a. According to theory (z 2 ) = 
1/2 for a — > and (z 2 ) = 1 when a — > 1. We compute 
(z 2 ) using Brownian simulations as in Table [J (solid curve) 
and with the QTM. For a = 0.8 and a — 0.9 the finite time 
simulations of the QTM did not converge, as shown in Fig. [5] 
However extrapolating the data (see Fig. [6| we get excellent 
agreement between simulations of the trap model and our 
theory based on weak sub-ordination breaking. 



where b\ and 62 are found in Appendix D. In the limit 
we find 



limg Q (0 : 
«->■! v2tt 



the expected Gaussian behavior. 



1 _e 

e 2 



(59) 



IX. CRITICAL SLOW DOWN a -» 1 

As discussed in [Z(| close to the critical point a = 1 
convergence of direct simulations of the QTM is ex- 
tremely slow. In contrast simulations of Brownian tra- 
jectories, using weak sub-ordination scheme converges in 
reasonable time, at least on our computer. In this sense 
our approach is much more efficient compared with di- 
rect simulation of the QTM. In fact we believe that our 
scheme is the only numerical tool available today for the 
investigation of the limit a 1. 

In Fig. [5] we show (z 2 ) versus a. (z 2 ) was obtained 
by several means: i) simulation of the QTM which give 
(x 2 ), and then with Eq. (l33l) we extract (z 2 ) (ii) Brown- 
ian simulations on a lattice (results in Table [J), and (iii) 
analytical theory Eqs. (j^2l IS"Tj) . For a > 0.8 our simula- 
tions of the QTM did not converge even for t = 10 9 . To 
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Subordination t=10 
□ Subordination t=I0 6 
Subordination t=[0 ? 

Large ^ 

Small 1 
— Analytical Prediction 




FIG. 6: Here we show that A defined in Eq. (|60p approaches 
zero extremely slowly when a = 0.8 and a = 0.9. Following 
the observation of Bertin and Bouchaud, close to the tran- 
sition from anomalous to normal types of diffusion we find 
critical slowing down. We see that A(t) ~ 0.1 at times 10 9 
when a = 0.8 and hence simulations did not converge at a 
time scale which is on the verge of our numerical capabilities 
(for a = 0.9 the situations is worse). By extrapolation the fig- 
ure shows that the asymptotic theory is reached albeit slowly. 
The dashed line is a guide to the eye with a t~ 01 behavior. 



FIG. 7: For the case a = 0.9 we present the scaling PDF 
obtained using weak subordination scheme Sec. |V] to- 



gether with analytical curve Eq. (|52|l and large and small £ 
expansions Eqs. (|55l [58}. Similar direct simulation of the 
quenched trap model do not converge even when t = 10 9 (see 
Figs. I5I6[) The figure demonstrates that with our approach 
based on weak sub-ordination breaking reasonable converges 
is found already for t — 10 5 . Hence close to the critical point 
a — 1 weak subordination breaking can deal with the critical 
slow down. 



check this issue better we define the deviation 



X. DISCUSSION 



(x 2 )aT 



A(i) = 



6&) 



1+c 



-)(t/AV°)" 



(60) 



According to Eq. ([33} Iim t ->.oo A(i) =0. In Fig. ©we 
present A(i) versus time. Here (a; 2 ) is obtained from 
QTM simulations and (z 2 ) from Brownian trajectories 
(see Table 1}. In Fig. M we show that A(i) ~ i Q_1 
for a — 0.9 and observe a very slow decay towards the 
asymptotic value A(t) — > 0. Simulations of the QTM did 
not converge (even for t — 10 9 ) however extrapolating 
the data (assuming nothing dramatic happens for times 
larger than 10 9 ) we can conclude that A(t) — > and in 
that sense our theory is consistent with the simulations. 

To over come critical slow down in the region a — > 1 we 
use weak subordination scheme Sec. [V] instead of direct 
numerical simulations of the QTM. In Fig. [7] we show 
the diffusion front. Good agreement between analytical 
predictions (valid for a — >• 1 ) Eqs. ([521 [551 155} and 
weak subordination breaking algorithm is presented for 
a = 0.9. 



The main focus of this paper was on the diffusion front 
(P(x,t)) of random walkers in the quenched trap model 
in one dimension. We showed that (P(x, t)) is found with 
a Levy time transformation acting on Brownian motion 
stopped at the operational time S a . Thus we map the 
random walk in disordered environment to a Brownian 
motion which is stopped at the a coverage time S a . This 
new type of Brownian motion is interesting on its own 
right. For example we have found a transition from a 
V shape to a Gaussian behavior for the scaling function 
B a (z) describing this motion. Properties of this function 
determine the statistics of diffusion in the QTM. For a 
close to 1 and we obtained analytical expressions for 
B a (z) and (P(x,t)} while numerical information easily 
obtained from Brownian simulations provide a detailed 
description of the diffusion front in the range < a < 1. 

For a — > our formulas reduce to the renormalization 
group results obtained by C. Monthus (HJ. The approach 
presented here is an alternative to the renormalization 
group method. Its advantage is that it is capable of deal- 
ing with the whole spectrum of a, at least numerically, 
including in the critically slowed down regime of a — > 1. 

Is our method general or is it limited to the one dimen- 
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sional quenched trap model? Clearly our approach can 
be extended to higher dimensions, or for random walks 
with biases. As mentioned in the introduction beyond 
the critical dimension, the QTM belongs to the domain 
of attraction of the CTRW. Hence for an ordinary ran- 
dom walk on a lattice in dimension three we expect that 
S a is non-random and equal to c a s when time s is large, 
and c a is a constant so far not determined. In that case 
usual subordination method works for the corresponding 
trap model. Hence once the constant c a and the dif- 
fusion coefficient of the corresponding discrete time ran- 
dom walk arc determined we have the statistical informa- 
tion needed for the determination the diffusion front of 
the QTM. Detailed analysis of the QTM for dimensions 
higher than one and for biased processes, using methods 
developed here, are left for future work. 



XI. APPENDIX A 

Here we summarize some known results on one sided 
Levy stable random variables, which are used all along 
this work. By definition i a) i i(t) is the inverse Laplace 
transform of exp(— u a ). The large t series expansion 

= i - r (i + H x s . n t _ (aB+1J< 

71=1 

_ (61) 
The asymptotic small t behavior is |43| 

l a ,i,i(t) ~ Bt-°e~ Kt ~ T (62) 

where 

1 — a 2(1 — a) 
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B = {[2tt(1 - a) 



1/2 



(63) 



Closed form PDFs are found by summing the series 
Eq. (pTTj) for specific choices of a |35l . |45| . For example 
we insert the series Eq. (|6"Tj) in Mathematica and use the 
command Simplify to get Levy PDFs in terms of a com- 
bination of Hypergeometric functions (e.g. for a = 1/4). 
Similarly Levy's PDFS with a = 1/4,1/3,1/2,2/3,9/10 
can be expressed in terms of special functions. In this 
way we construct stable distributions. Some care must 
be practiced since in some occasions we found that for 
extremely small t Mathematica yields wrong results [35j . 
This problem is easily fixed since we can use Eq. (|62[) 
in that regime. Further the problem is not crucial in the 
sense that it is found for so small t that practically the 
PDF there is zero, though if one is not aware of this issue 
solving integral transformations like Eq. (|4"9"]l can lead to 
wrong results. A useful special case is a = 1/2 since 



(64) 



Chambers et al [59| show how to generate a stable 
random variable we call r\ from a one sided Levy PDF 



sin((l- a )fl)(sinM)r /(1 - a) 

a ( & ) : . „,i/(i-a) ' < < 7r - 



(sin 6>) 



(65) 



Then 77 = [a(9)/W]^ 1 ^ a ^ a where 9 is a uniform random 
number on (0, it) and W is a random variable drawn from 
a standard exponential distribution: W — — ln(x) were 
x is uniform on (0, 1). 



XII. APPENDIX B 



In this Appendix we obtain the distribution of 77 de- 
fined in Eq. (|15[) . We are interested in random walks with 
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fixed S a and in the limit S a — > oo. A large S a implies 
also a large number of steps (denoted with s), however 
since S a is fixed s remains random. Our starting point 
is Eq. (HU) 



ing time PDF is in the small u limit 



i>(u) = 1 - Au a + Bu p + 



(67) 



( e -™'} = n~ 



n x u 



(66) 



It is important to notice that the visitation numbers {n x } 
are determined by the probabilities of jumping left and 
right (equal 1/2 in our model) and that these random 
numbers do not depend on the waiting times since here 
S a is fixed. Hence statistics of these visitation numbers 
are determined by simple binomial random walks. 

The Laplace r — ¥ u transform of a rather general wait- 



where as mentioned 0<a<l,A>0 and (3 > a. The 
goal is to show that when S a — > oo parameters like B 
and j3 are not important. To see this insert Eq. (|67p in 
Eq. flBBJ and find 



IT 



B v ; v" 



This can be rewritten as 



(68) 



OO OO 



A 2 (n x ) a (n y ) a 2a , u Y,Z-oo(n*y 



-^) = i-Au a + y y — VLx > y 'r u 2a +B ^=- °°; x */+•••. m 

^ 2 (s r k 



We note that 



oo oo oo oo 



hence 



E E Kr(%r= e Ew>/"Ew 2a = ( 5 «) 2 -^ (™) 

x— — oo y——oo,y^x x — — oo y— — oo x— — oo 



(e )-! Aw + 2 (5q)2 w + jB ( 5q )/V« u + ■ (71) 



We use S 2a /(S a ) 2 — > which is justified at the end of 
this Appendix and hence 

{S a f - S 2a . . 

when S a -> oo. Similarly Sp/(S a )P/ a -> for a < /3. 
Summarizing we find 

(e-" u ) ~ 1 - Au a + + • • • = e- Au ° (73) 

The parameters B and /3 are unimportant. Further for 
a typical path there is no trace of the random variables 
{n x } in the final expression Eq. (|T3|) . The latter Eq. 
implies that the PDF of 77 is a one side Levy PDF as 
stated in Eq. (fT7|) . 

To better estimate the convergence to this law we use 
a result obtained in Appendix C. There we show that for 
a binomial random walk with s steps we have 

(S a ) = C a s {1+a ^ 2 (74) 

where C a is a constant Eq. ([§(T]) . We then assume the 
following relation to hold 

S a = rs ( 1+ «)/ 2 (75) 



where r is a random variable which is independent of 
the number of steps s. Since in the QTM jumps to 
nearest neighbors have probability 1/2 like the bino- 
mial random walk and since the statistics of the vis- 
itations numbers {n x } are independent of the waiting 
times (for fixed S a ) we may use Eq. (|74p derived for 
the binomial random walk to analyze the QTM. We have 
S 2a oc s ( 1 + 2q )/2 and hence S 2a <x (S a )( l + 2a V(. 1 +°>) so 
S 2a /(S a ) 2 oc (S a )~ 1 ^ 1+a * > which goes to zero in the scal- 
ing limit S a ~ > 00 as we stated . Similarly Sp / {S a )^^ a oc 

s (l+0)/2/ s /9(l+a)/2a = s (a-/3)/(2a) which ap p r0 acheS zero 

since a < f3. 



XIII. APPENDIX C 

We consider a binomial random walk on a one dimen- 
sional lattice. Time s is discrete s = 0, 1,2, ... and the 
particle has probability one half to jump to its nearest 
neighbors on its left or right. The walk starts on the ori- 
gin x = 0. We now calculate the average (S a ) for an s 
step random walk. For that aim we obtain ((n x ) a ) and 
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then sum over x 



(S a ) = £ <(n x )«) 



(76) 



X — — oo 



We consider this problem in the continuum limit of the 
model, namely we consider Brownian motion (see details 
below). Thus our final expression for (S a ) describes the 
limit of large s. 

Let P s ,x{n x ) be the probability of making n x visits on 
lattice point x in the time interval (0,s). As usual in 
these problems the Laplace transform 



Pu, X (^a:) 



e- us P SiX (n x )ds. 



(77) 



is useful. Here we already started taking the continuum 
limit, since in the discrete time random walk s is not a 
continuous variable. We avoid a formal transition from 
a discrete random walk to the continuum limit, to save 
space and time. 

For a random walk starting on the origin let r be the 
first time the particle reaches lattice point x and / x (t) 
its PDF. r is a first passage time for an unbiased ran- 
dom walk and it distribution is well known [311 ] . From 
symmetry /_ x (t) = / x (t). The number of visits on lat- 
tice point x, n x , is determined by a first passage time 
from the origin to point x, and then by the probability 
to revisit point x. Due to translation symmetry of the 
random walk the probability of n x — 1 revisits to a lattice 
point x, in a time interval s — t (once reaching that point 
at r) is identical to the probability of n x — 1 visits on 
the origin (starting on the origin) within the same time 
interval. Hence translation symmetry gives 



Ps,x \P*x ) 



f x (T)P s - Tfi (n x - l)dr. 



(78) 



Using convolution theorem P u , x {n x ) — f x (u)P U fl(n x — l). 
Here P s ,o(no) is the probability to visit the origin no 
times in the time interval (0, s). 
For the origin x = we have 



Pu, x =o( n o) 



LAM 



h(u) 



ii-AW 



n = 
n ^0. 



(79) 



Here fi(u) is the Laplace transform of /i(t). To reason 
for Eq. (|79p note that if no = we have P s . x =o( n o = 0) = 
1 — J /i(r)dr which is the probability of not returning 
to the origin. To see this note that after one jump the 
particle is either on i = 1 or i = - 1 and hence for 
no to remain zero the particle must not return to origin 
(of-course /i(r) is the PDF of first passage times from 
x = 1 or x = — 1 to the origin). Applying the convolution 
theorem theorem of Laplace transform to P s , x =o{no — 
0) = 1 - J s /i(T)dr we get the first line in Eq. (1791 . 
Note that the original stay on the origin; at time s = 0, 
is not counted so we may have no = once the particle 



never returns to origin. The probability that no = 1 is 
given by 

P.,i(no = l)= / /i(r)P s - T ,o(no = 0)dr. (80) 
Jo 



Again using the convolution theorem we find Eq. (|79[) 
for no = 1. Similarly for n > 1. Using Eq. ([75]) it is 
easily shown that for x =/= 

Pu,M = { . - u n /™ ( 81 ) 



The Laplace transform of the first passage time PDF 



is 



f x (u) =exp(-V2^ 1/2 ) 



(82) 



The y/2 comes from the fact that the diffusion constant 
is equal 1/2, since the variance of jump lengths is unity. 
In time t, f x (r) is the one sided Levy PDF with index 
1/2 [see Eq. §4§] and / x (r) - r" 3 / 2 as well known [3lj]. 

We now calculate ((no)") u : the Laplace transform of 
{( n o) a )s 



((n Q ) a )u = / (n ) a P u ,o(no)dn (83) 
Jo 

where the integration (instead of summation) implies 
that we are considering the continuum limit of the ran- 
dom walk (i.e. Brownian motion). Inserting in Eq. 
Eqs. ((79152")) we find 



(84) 



((no)°)u~ (-7=) r(l + a) 



which is valid for small u corresponding to large s. Us- 
ing Eqs. (|81I82[) we get the a moment of the visitation 
number for lattice point x 

(K) a ) u = e- lA1/2 r(l + a)fij u- 1 -^ 2 . (85) 

Notice that ((n x ) a ) u ~ f x (u)((no) a ) u reflecting the first 
arrival at x and the translation symmetry of the lattice. 

Denote (S a ) u as the Laplace s — > u transform of (S a } s . 
In the Brownian limit we replace the summation in Eq. 
(|76p with integration 



(S Q ) U = 2 / ((n x ) a ) u dx, 
Jo 



(86) 



and with the help of Eq. ([8"5]l 

(5 a ) u ~V2 1_a r(l + a)u- a / 2 -«/ a . (87) 
Using the Laplace pair 



S 2 
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we find with simple identities for the Gamma function 
[6l| , the main result of this Appendix 



(S a ) = C a s i° 



with 



r 



2 2 * 2 r(i 



(89) 



(90) 



For a = 1 we have — y~!°° n r = s a result that is 
retrieved from Eq. (|89p since C\ — 1. In the limit a = 
we have Sq equal to the span of the random walk, namely 
to the number of distinct sites visited by the walker. As 
mentioned in the main text, in this limit we retrieve Eq. 
(|2"61 found a long time ago in (60| . 



XIV. APPENDIX D 

Here we investigate the function g a (£) in the limit 
where a < 1 is close to unity with the Gaussian ap- 
proximation for B a (z). Changing variables in Eq. (|49j) 
according to S a = t a rj~( 1+a ^ 2 we find using the defini- 
tion in Eq. @ 

g a (0 = 

(91) 



V(l+") 



2a 



where 



e (c a ) 



2/(l+a) 



/2- 



(92) 



Eq. (|9"Tj) is a Laplace transform. Inserting in Eq. (f9"T|) 

£ = 0, changing variables according to y = ^(i+ a )/( 2a ) 
and using Eq. flU) we get ffct (£ = 0) Eq. (J5BJ(. The 
small m limit (corresponding to small £) of Eq. (j9"Tj) is 
controlled by the large rj behavior of 



-l-a/2 



(93) 

where we used the large 77 expansion of stable PDFs Eq. 
(ItJTI) . Using the Tauberian theorem, noting that it"/ 2 
and t ? -( 1+q / 2 ) /|r(-a/2)| are Laplace pairs, Eqs. (|91I93I) 
together with some identities of the Gamma function [61] 
give Eq. (|55|) . 

In the opposite limit of large £ (i.e. large u) we use 
the small rj behavior of ^ 1/(2q) Z q ,i,i ( ?? ( 1 + Q )/( 2tt )) in Eq. 
(|9T|). For that aim we use the small ij behavior of one 
sided stable PDFs Eq. We get 



g a (0 = C 



dr/ 



(94) 



where k is defined in Appendix B Eq. (|63[) . 

3 — a 1 + a 

7 = 4(1 - a)' = 2(1 - a) 



C = B 



{C a ) 1 « 1+a l + a 
v/2^ 2a 



(95) 



and B is defined in Eq. (|63|) . 

We now use steepest descent method. Let ^1(77) = 1*77 + 
nrj^ s . The extremum is on rj e which is determined us- 
usual from dh/drj = so r\ e = (u/k5) _1 /' 1+<5 ''. Using the 
expansion 



-5-2a2 



(96) 



where A = 7/ — r] e is small. We then have after extending 
the domain of integration 

g a (0 ~ 



C (?7 e )~ 7 e~ hM / exp[-5(l + 5) K (r 1e )- s - 2 A 2 /2]dA. 

J — OO 

(97) 

Solving this Gaussian integral 

g a (0 = du-»e--™ 5ni+5) , (98) 

where C = CV2^[5{\ + S)k]-^ 2 (kSY 1 , k = (k < 5) 1 /( 1 +' 5 ) + 
k(k5)- 5 /( 5+1 ' and fj, = (1 - a)/(3 - a). Changing back 
to the variable £ (instead of u) using Eq. (|92p we get Eq. 
(El 



C 





9 


»(0 = 








J 






-1 


2 


exp < 




2 





(99) 



To prepare for the limit a — > 1 we rewrite 



C = (C a )~ 
Using lim a _j.i C a = 1, 
lim 



i_ 1 + a 



2a ^ (1 + 5) K 



(100) 



/2vr 



(101) 



and lim Q _j.i ft = 1 we find the expected Gaussian behav- 
ior Eq. (|59|) . Rewriting 



(102) 



when £ >> 1 with 61 = ^(1 + a)/[27ra(3 - a)]D, b 2 
[(3 - 2a)/2] # 2 and £>=[(! + o) 1 - a a a C a ] 1/(3 " a) . 
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